Electronic and magnetic properties of the ionic Hubbard model on the striped 

triangular lattice at 3/4 filling 
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We report a detailed study of a model Hamiltonian which exhibits a rich interplay of geometrical 
spin frustration, strong electronic correlations, and charge ordering. The character of the insulating 
phase depends on the magnitude of A/\t\ and on the sign of t. We find a Mott insulator for 
A ^ U ^ \t\; a. charge transfer insulator for L*' ^ A ^ \t\; and a correlated covalent insulator for 
U ^ A ~ \t\. The charge transfer insulating state is investigated using a strong coupling expansion. 
The frustration of the triangular lattice can lead to antiferromagnetism or ferromagnetism depending 
on the sign of the hopping matrix element, t. We identify the "ring" exchange process around a 
triangular plaquette which determines the sign of the magnetic interactions. Exact diagonalization 
calculations are performed on the model for a wide range of parameters and compared to the strong 
coupling expansion. The regime [/ ^ A ~ |f| and t < is relevant to Nao.5Co02. The calculated 
optical conductivity and the spectral density are discussed in the light of recent experiments on 
Nao.5Co02. 

PACS numbers: 71.10.Fd, 71.15.-m,71.27.+a 
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I. INTRODUCTION 

Many strongly correlated electron materials exhibit 
a subtle competition between different magnetic and 
charge ordered states, and between metallic, insulating, 
and superconducting phases. Widely studied (and poorly 
understood) materials include cupratc superconductors fi 
organic charge transfer salts,— manganites with colossal 
magnetoresistance?^ heavy fermion compounds;^ and the 
iron pnicitide superconductors A fundamental theoret- 
ical challenge is explaining the hierarchy of energy scales 
and competing phases in these materials. The energy 
scales (such as the bandwidth and Coulomb repulsion) 
associated with the relevant electronic orbitals (and mi- 
croscopic Hamiltonians such as Hubbard models) arc typ- 
ically of the order of eV. In contrast, the energy scales as- 
sociated with the temperature and magnetic field depen- 
dences of transport properties and energy differences be- 
tween competing phases are often several orders of mag- 
nitude smaller. Frustration of spin or charge ordering by 
competing interactions due to the geometry of the crys- 
tal lattice can enhance these effects. In addition, it is 
not clear what physical changes arc produced by chemi- 
cal doping. For example, does adding charge carriers just 
change the band filling or are there significant effects due 
to the associated disorder and changes in the electronic 
structure? 

Here we report a detailed study of a specific strongly 
correlated electron model, the ionic Hubbard model on 
the triangular lattice at 3/4 filling with a stripe potential. 
The model illustrates how the interplay of geometric frus- 
tration and strong correlations lead to competition be- 
tween different magnetic orders, charge ordering, metal- 
lic, and insulating behaviours. One concrete realisation 



of the model is that it may be the simplest many-body 
Hamiltonian that can describe Nao.5Co02.— "^"^ Elsewhere 
we have reviewed experimental results on this material 
and described recent theoretical attempts to describe its 
unusual propertiesi^ When the filling x in Naa;Co02 is 
close to other commensurate values, such as 1/3,2/3, or 
3/4, the system is still described by an ionic Hubbard 
model but the on-site potential has a different form and 
commensurability, depending on the ordering arrange- 
ment of the sodium ions,^ At incommensurate values of 
X one expects phase coexistence of multiple Na-ordering 
phases 1^ 

The rest of the paper is organized as follows. In Sec. 
im we introduce an ionic Hubbard model on a triangu- 
lar lattice including a discussion of its phase diagram. 
In Sec. mil we analyse the model's ground state prop- 
erties using the Lanczos exact diagonalization technique 
on finite-size clusters. Dynamical properties such as the 
spectral density and optical conductivity are discussed in 
Sec. IIVI Finally, a summary of the main results and their 
relevance to Nao.5Co02 is given in Sec. |Vl We have also 
studied the same model using a complementary method, 
mean-field slave bosonsi^ At appropriate places in the 
paper we compare and contrast the results. 



II. THE IONIC HUBBARD MODEL ON A 
TRIANGULAR LATTICE 



The Hamiltonian of the ionic Hubbard model is 

H = -t (c+Cj^-hc^;c,^) + C/^n,|n,|-f ^e, 



(1) 
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FIG. 1: (Color online) Spin and charge order in the 3/4-filled 
Ionic Hubbard model ([TJ on a triangular lattice in the limit 
U » A » \t\. A and B denote the inequivalent sites of the 
lattice. C-type antiferromagnetism (left) is found for t > 
in contrast to G-type antferromagnetism (right) for t < 0. 
The exchange couplings J and J± are defined in Eq. ((Tjl for 
the appropriate t — J-model in Eq. (|6]). A ferromagnetic 
exchange coupling, J, between neighboring A-sites occurs for 
the parameter range: < 5t < A < V2U. 



where creates an electron with spin a at site i,t is the 
hopping amplitude between neighbouring sites and U is 
the effective on-site Coulomb repulsion energy between 
two electrons. We set = A/2 for the A-sites and e,; = 
—A/2 for B-sites (cf. Fig. [T]). The A-sites form rows 
which alternate with the B-sitcs of the triangular lattice. 
By a particle-hole transformation: hi„ model ([T]) 

becomes a 1/4-fillcd (with holes) ionic Hubbard model 
with the sign transformation: t —t and A —A. 

In what follows we will discuss the charge gap which is 
defined for the model on a finite lattice with N electrons 
and 2A^/3 lattice sites by 

A, = Eo{N+l) + Eo{N-l)-2Eo{N), (2) 

where Eo{M) is the ground state energy of the system 
with M electrons. 

A. Previous theories of the ionic Hubbard model 

The ionic Hubbard model ([T]) on the striped triangu- 
lar lattice of Fig. [1] considered here contains geometrical 
frustration. This is in contrast to most previous work 
which has focused on bipartite, i.e., unfrustrated, lattices 
with different site energies, A, on each of the bipartite 
sublattices and at half-filling. Some of the interest in this 
model can be appreciated from the half-filled atomic limit 
{t = 0). For U > A the charge gap for the addition of 
particles defined in Eq. ([2]), Ac = U— A; thus the system 
is a Mott insulator. But, for U < A, Ac = A — U ; and the 
system is a band insulator. At the point U = A this gap 
vanishes. Therefore a key question is what happens at the 
band to Mott insulator transition away from the atomic 
limit; in particular what happens to the gapless point - 
does it expand give a metallic phase? Further impetus 
comes from the proposals that the ionic Hubbard model is 



important for understanding ferroelectric perovskites^ 
organic charge transfer salts ^iii^ transition metal oxide 
heterostructureSfi^ and non-linear electronic polarizabil- 
ity in transition metal oxides and, as discussed below 
the rich electronic phases observed in Aa;Co02 [A=Na, 
K, Rb] (Ref. i). 

The most studied case is the half filled onc-dimensional 
chain with different site energies for odd and even num- 
bered sites. This model shows three distinct insulat- 
ing phases: a band insulator; a (ferroelectric) bond or- 
der wave insulator; and a Mott insulator.— Metallic be- 
haviour appears to be limited to the point in the phase di- 
agram where band insulator gives way to the bond order 
wave insulatorii^ Continuum limit bosonization calcula- 
tions suggest that adding a next nearest neighbour hop- 
ping, t', (which is equivalent to studying the zigzag chain) 
induces a large metallic region in the phase diagram^ 
suggesting that even in one-dimension frustration already 
plays an important role in the ionic Hubbard model. 

The infinite-dimensional ionic Hubbard model 
has been studied using dynamical mean-field theory 
(DMFT)^ii^ which treats the on-site quantum dynam- 
ics exactly but ignores spatial correlations such as those 
associated with antiferromagnetic exchange. These 
papers studied the, bipartite, Bethe lattice and found 
that a metallic phase separates the band insulating 
phase from the Mott insulating phases, in, at least, some 
parts of the phase diagram. 

In two-dimensions most previous work has focused on 
the half-filled square lattice with site energies alternating 
in a checkerboard pattern. This model has been studied 
using both cluster DMFT— and determinant quantum 
Monte Carlo (DQMC).^^'^ These studies all suggested 
that a phase with non-zero spectral weight at the Fermi 
energy exists between the band insulator and Mott in- 
sulator phases, at least in some of the phase diagram. 
However, the has been some debate over whether this 
phase is metalli o^°'^^ or bond orderedJ^ 

There has been far less work on the ionic Hubbard 
model away from half-filling. However, Penc et al^ stud- 
ied the quarter-filled ionic Hubbard model on the zigzag 
ladder and found a competition between ferromagnetism 
and a paramagnetic phase with strong antiferromagnetic 
correlations. Bouadim et ai^^ studied the ionic Hub- 
bard model on a square lattice with a checkerboard po- 
tential across all possible fillings with DQMC. The most 
interesting features they found, away from half filling, 
were Mott insulation at quarter and three quarters fill- 
ing (which are related by the particle-hole symmetry of 
this model). Bouadim et al. did not find any evidence of 
magnetic order in this phase. However, as the magnetic 
interactions arc 0{t'^/UA'^) the absence of magnetic or- 
der may be due to non-zero value of the temperatures 
they studied. 

As well as bipartite arrangements of the different site 
energies there has also been considerable interest in ran- 
dom arrangements of site energies. Laad et al^ stud- 
ied a system with a gaussian density of states and a bi- 
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modal distribution of site energies in infinite dimensions 
for various impurity concentrations, n, and fillings, I — S. 
They found insulating states for S = 1 — n and suffi- 
ciently large A and U. Byczuk et al. have used DMFT 
to study the (frustrated) fco2^ and (bipartite) Bethe^i 
lattices in infinite dimensions with a bimodal distribu- 
tion of site energies and half the sites taking each value 
of the site energy. At one-quarter filling they found that 
metal insulator transitions occur on both lattices when 
both U and A are sufficiently large. They also noted 
that DMFT does not capture some of the possible effects 
of the disorder, such as Anderson insulating phases. In 
two dimensions Paris et al^ used DQMC to study the 
square lattice at a range of fillings with 1/8 of the sites 
randomly chosen to have a different site energy than the 
rest of the lattice. They found that this model displayed 
Mott insulating, band insulating, Anderson insulating, 
and metallic phases. 

Marianetti and Kotliar— simplified our suggestion^ of 
that Eq. ID) is the appropriate effective low Hamiltonian 
for Naa;Co02 by further assuming that Na-ordering is 
of secondary importance and hence treated the potential 
due to the Na ions as random. They then used density 
functional theory to show that the distribution of Co site 
energies is bimodal and to parameterise the Hamilton- 
ain ^ ioi X ~ 0.3 and 0.7. Finally they calculated the 
high temperature (T > 100 K) susceptibility for these 
dopings and found them to be in qualitative agreement 
with experiment. We deal here with the case x ~ 0.5 for 
which Na-ordering of the stripe- type (cf. Fig[T]) has been 
observed in experimentsi^i^ 

B. The non-interacting model {U — 0) 

For U = model ([1]) can be diagonalized straightfor- 
wardly leading to two bands, denoted ±. We introduce 
creation and destruction operators: 

(i±a = "kT(cika + ^kT4ka)' (3) 

where c^j^^ and c^j,^ act on the Bloch states associated 
with the A and B sublattices, respectively, and 

A/2± y^A2/4-h (4icos(fc^/2)cos(fc.yV3/2))2 
At cos(fc^/2)cos(fc.yV3/2) 

{±) 

with the normalization constant: ak± = 1/ -^/l + |Ak±p- 
The energy dispersion of the two bands is: 



e±(k) = -2tcosfcj;±yA2/4+ (4t cos(fc^/2) cos(fcy V3/2))2, 

(5) 

with kx and ky defined in the reduced (1 x VS) Brillouin 
zone with lattice parameter a = 1. At 3/4- filling and 
for any A, there is always at least one band crossing the 
Fermi energy and so the system is metallic. The -I- band 
is half-filled and the — band filled for t > and A > 0, 
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FIG. 2: (Color online) Schematic phase diagram of the Ionic 
Hubbard model[l]on a striped triangular lattice at 3/4-filhng. 
The transition lines are based on the lowest order corrections 
in a strong coupling analysis and on Lanczos diagonalization 
calculations. The t < (top panel) case is relevant to the 
Nao.5Co02 insulator. Insulating phases at strong coupling, 
U > W , of different types are found ranging from a charge 
transfer insulator (CTI), a Mott Insulator (MI) and a cova- 
lent insulator (CI). The bandwidth of the model for A = 
(the isotropic triangular lattice), W = Q\t\, is effectively re- 
duced to 4\t\ corresponding to one-dimensional chains as A 
increases. At exactly A = oo, the system is insulating for any 
nonzero U as expected for a half-filled Hubbard chain due to 
Umklapp processes. The blue line is an estimate of the crit- 
ical U for the metal-to-insulator transition which follows the 
effective bandwidth dependence with A. Insulating phases 
for t < display G-type antiferromagnetic (AFM) (see Fig. 
[T]) correlations whereas a C-AFM region (see Fig. [T]) for t > Q 
is obtained from the condition, J < 0, to Equation ((Tjl which 
we assume valid for U > 9|f|. The marked A = axis for 
t > Q and above U « 5|f| indicates the occurrence of ferro- 
magnetism as predicted by DMFT of the Hubbard model on 
a isotropic triangular lattice^. 



whereas for t < Q this only occurs for A > 0.64|t|. For 
A there is only one band, which has a width of 
W = 9|i|. 
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C. Phase diagram 

The ionic Hubbard model ([T]) on a triangular lattice 
contains a rich phase diagram resulting from the inter- 
play between geometrical frustration, strong Coulomb re- 
pulsion and charge ordering phenomena. A schematic 
phase diagram can be constructed by first considering 
some simple limits: 

(i) t = (atomic limit): all A sites are singly occupied 
while B sites are doubly occupied. The system is insulat- 
ing with a charge gap: Ac = min(A, U). For A > [/ it is 
a Mott insulator (MI) with Ac = U, whereas for A < C/ 
it is a charge transfer insulator (CTI) with Ac = A. 

(ii) U = (non-interacting limit): As discussed above 
the model is always metallic regardless the value of A. 

(Hi) A = 0: For any U the model reduces to the regu- 
lar Hubbard model on the isotropic triangular lattice at 
3/4 filling. For large C//|i|, it is equivalent to the t — J 
model on the triangular lattice. Dynamical mean-field 
theory calculations^!^ give a ground state that is metal- 
lic. DMFT— calculations on the Hubbard model for large 
U and variational Monte Carlo calculations on the t — J 
model give this metallic ground state as paramagnetic 
(ferromagnetic) for t < (t > 0)3 

(iv) A = oo: As the B-sites can be completely pro- 
jected out from the Hilbert space, the model is mapped 
onto decoupled half-filled Hubbard chains. Hence, the 
system is (Mott) insulating for any non-zero positive U; 
i.e., there is a charge gap Ac ^ 0, and there are antifer- 
romagnetic correlations (with power law decay) and no 
spin gap. 

(v) U » A >> |<| 7^ 0; For finite but small t, virtual 
hopping processes lead to effective magnetic exchange 
couplings between the A sites. The effective low-energy 
t — J — Jdiag — J± Hamiltonian for the holes is 



H 



t J2 PiKahj. + h+MP + J ^ [Si . Sj 

Jdiag ^ I^Si • Sj 



(y) 



(6) 



where {...}, and [...] denote sums over intra- A-chain, 
inter- A-chain sites, respectively. The sum over (...) is be- 
tween an A and nearest-neighbour B-sites. The projector 
P = Hi [1 — Ui^Uii] forbids double occupation of holes on 
any lattice site. The dynamics of the electron-doped sys- 
tem relevant to Naj;Co02 with electron occupation l-\~x 
is related to the hole-doped system, with filling 1 — x, 
through the replacement: t —t and — > — e.^ leaving 
the exchange parameters unchanged. 

The exchange couplings J, J± and Jdiag can be 
obtained through a strong coupling expansion using 
Ralcygh-Schrodingcr perturbation theory on the hopping 
term around the configuration in which all B sites are 




FIG. 3: (Color online) Cluster shapes of different sizes used 
in exact diagonalization calculations. 



doubly occupied and A-sites singly occupied (see Ap- 
pendix 1^ . This leads to an effective exchange coupling 
between electrons in A-sites in the horizontal direction: 



_ At^ 8t^ 



16t^ 
AC7 



Oit^ 



(7) 



and in the perpendicular direction 



J± 



A2 



1 

U 



1 



2A + U 



1 

2A 



+ 0(i^). (8) 



The exchange coupling Jdiag between A and B sites is 
1 1 



Jd 



lag 



2t^ 



U+A U-A 



+ 0(<^), (9) 



which is blocked if the B sites are doubly occupied but 
recovers the correct At'^/U exchange interaction as A ^ 
0. 

The second and third terms in J arc antiferromagnetic 
(AFM) for t < and ferromagnetic (FM) for t > 0. 
Higher order contributions to J are AF and can be found 
in Appendix [Xl 

A schematic phase diagram of the 3 / 4-filled ionic Hub- 
bard model on a triangular lattice ([1]) is shown in Fig. 
m The transition lines are extracted from the limits (i)- 
(v) discussed above and exact diagonalization calcula- 
tions for intermediate parameter regimes. Apart from the 
Mott insulator (MI) and charge transfer insulator (CTI), 
our numerical analysis suggests the presence of a cova- 
lent insulator (CI) in the range A ~ 0(|i|) and U » \t\. 
Depending on the sign of t, different spin arrangements 
occur as shown in Fig. [TJ The condition J = separates 
AF from the FM region which occurs in the parameter 
range : 5t < A < Vm and is plotted in Fig. H 



D. The model on two and four-site clusters 

In this section we explore the nature of the ground 
state of the model ([l} on two and four-site clusters. The 
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FIG. 4: (Color online) Ionic Hubbard model on two- and four- 
site clusters. The energy level diagram for 3 electrons on two 
sites (left) and the four-site cluster with six electrons (right). 



two-site cluster incorporates charge ordering phenomena 
in the presence of on-site Coulomb interaction. The four- 
site toy model also contains geometrical frustration ef- 
fects present in the full model Ground state prop- 
erties of the clusters are discussed in terms of valence- 
bond (VB) theorjs^ when appropriate. For t < and 
U >> |t| the ground state wavefunction is accurately de- 
scribed by the resonance between possible valence bonds. 
Our analysis indicates that the charge gap of the clusters 
is enhanced with A due to the differences between the 
two-electron and three-electron bonds between the dif- 
ferent A and B sites. 



1. Two-site cluster 

We first consider three electrons in two incquivalcnt 
sites (one A, the other B) separated by an energy A. 
The energy levels for this cluster are sketched in Fig. H) 
The Hamiltonian is: 

+ A/2{nA~nB). (10) 

For U = 0, the charge gap of the cluster is Ac = 
for any A due to the degeneracy of the ground state. In 
order to evaluate the dependence of the gap on A we 
first obtain the ground state energies with N =2, 3 and 
4 electrons: 



So (3) 
So (4) 



-2r 

2U. 



1 



1 



C/-f- 



U-A 



,u» \t\ 



■4*2 



(11) 



In the limit A ^ and U >> |*|, the charge gap of the 
cluster is 



Ac«2|i| 



At' 



AM - - ^ - ^ 



At 



At' 



The first contribution to Ac is present even for A = as 
expected from the bonding-antibonding splitting of the 
cluster and will go to zero in the infinite system. The 
term proportional to A^ comes from the different depen- 
dances of £'(3) and E{2) on A: E{2) has a weaker depen- 
dence than E{i). This is due to the different nature of 
the two-electron and the three-electron bond. The former 
is accurately described by a correlated VB between an 
electron on an A site and an electron on a B site whereas 
the latter is described by a single hole in an antibond- 
ing "molecular" orbital. The two-electron and three- 
electron bond energies (the energy needed to break a 
bond between inequivalent sites) are A/2 — VA^ -f 4*^/2 
and Eo{2), respectively. Hence, the two-electron bond 
becomes weaker with A while the three-electron bond 
is strengthened with A. This is known from quantum 
chemistry^ and can be understood as being a conse- 
quence of the presence or absence of Coulomb repulsion 
between electrons. 



2. Four-site cluster 

We consider the four-site cluster (Fig. 0]) with two A- 
sites shifted by -l-A/2 and two B-sites shifted by —A/ 2. 
The cluster shown contains TV = 6 electrons (correpond- 
ing to 3/4- filling). The model Hamiltonian in this case 
is 

+ + A/2 ^ [n^-rij). (13) 

i ieA,j£B 

We first discuss the A ^ limit. The exact ground 
state energies of the fully frustrated t' = t cluster for 
= 5, 6 and 7 electrons are plotted in Fig. [5] for 
U ~ 100|t|. Consistent with the results for the two-site 
cluster we find that, Eq{5), has the weakest dependence 
on A of all the groimd state energies. The ground state 
wavefunction for A^ = 6 can be well described in terms 
of resonating valence bonds as shown in the Appendix 
IbI In contrast, the wavefunction for N = 7 consists of a 
single hole hopping around the cluster and so contains no 
Coulomb interaction effects. In this case a "molecular" 
orbital with a single hole describes the cluster and its en- 
ergy, Eo{7), has the strongest dependence with A. The 
different behavior of Eo{N), Eo{N + 1) and Eo{N - 1) 
is responsible for the increase of the charge gap with A 
as shown in Fig. [51 This behavior is in contrast to the 
U = case plotted in the figure. 

Effects of frustration. We now discuss the four-site 
cluster with t' ~ 0. In this case analytical formulas may 
be obtained for U = oo. The sign of the hopping is 
irrelevant here in contrast to the t' ^ t case. The ground 
state energies are: 



U 



(12) 



So (6) 



6 
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-12 



1 1 1 1 1 


1 ' 1 ' 
t>0 


r — Eo(7) 




- — Eo(5) 




- — Eo(6) 

1 1 1 1 1 


1,1, 



10 



A/|t| 




FIG. 5: (Color online) Dependence of ground state energies 
with A for the four-site cluster with t' = t and U = 100|t|. 
The Coulomb interaction U, 2U and 3U have been sub- 
stracted from the total energies -Bo (5), Eo{6) and Eo{7), re- 
spectively, for convenience. 



FIG. 6: (Color online) Dependence of the charge gap on A 
for the four-site cluster. Note the increase of Ac for any A for 
U >> \t\ of Fig. ((5)) and also how the dependence for A > |t| 
is quite different with A < \t\. 



III. GROUND STATE PROPERTIES OF THE 
^4^2 + A2 IONIC HUBBARD MODEL ON A TRIANGULAR 

LATTICE 



So(7) = -t-^ ^ . (14) 



In the limit A 0, wc find Eo{7) w -2t - ^ and 

£:o(6) w "2^/2^-^^^. Thus, £;o(6) has a weaker depen- 
dence on A than Eo{7), similarly to the fully frustrated 
t' = t cluster. 

Our small cluster analysis indicates that the charge 
gap, Ac, increases with A due to the different natures 
of the two-electron and three-electron bonds formed be- 
tween inequivalent sites. This result is not affected by 
the presence of frustration in the cluster at the qualita- 
tive level. However, geometrical frustration (t ~ t') leads 
to qualitatively different magnetic properties for different 
signs of t in contrast to the unfrustrated {t' = 0) case. 



Intermediate parameter regimes are explored based on 
Lanczos diagonalization on finite clusters with Ns=12, 16 
and 18 sites and periodic boundary conditions. Different 
cluster shapes have been benchmarked against the exact 
solution of the non-interacting model ([T]) and are shown 
in Fig. [3l The vectors defining the clusters are: Ti = 
niiai-|-ni2a2 and T2 = ?T.2iai -l-ri22a2, where nu and n2i 
are integers. A straightforward finite size scaling analysis 
is not possible because of the complicated changes in the 
cluster shape as the lattice size increases. 

We present results of the dependence of the charge 
order parameter, the charge gap and the spin correlations 
on A. Numerical results are compared to the weak and 
strong coupling limits as appropriate. 
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A. Charge order 

The charge order parameter is first computed for U ~ 
and compared to exact tight-binding results on the infi- 
nite lattice. This serves to calibrate the importance of 
finite-size effects on a cluster. Second, the effect of U on 
charge ordering is analyzed in detail. 

The charge order parameter is 

nB-nA = ^(vE-oKc^k^Ciska - c\^„CAk„)\'i'o) , (15) 

k,(T 

where \^o) is the ground state of the Hamiltonian H]). 

In Fig. [7|the charge order parameter is plotted for both 
signs of t on the Ns = l8 site cluster of Fig. [3] and com- 
pared to the tight-binding result {U = 0) of the extended 
system. Since this cluster gives the best agreement with 
the infinite limit of the non-interacting model, among all 
of the clusters that we studied, we mostly show results 
for this cluster in the rest of this paper. 

The effect of Coulomb repulsion on charge transfer is 
also shown in Fig. [T] The qualitative dependence of 
charge transfer remains unchanged with U. However, 
increasing U does suppress ng — a little for small and 
moderate A. 



In Nao.5Co02 the strong Coulomb interaction and 
weak charge transfer implyli^: U >> A and A ^ \t\. 
Note also that in this parameter regime charge trans- 
fer between A and B sites is weak: ub — tla < 0.2, for 
A ~ |i| which implies that a charge transfer insulator 
formed by doubly occupied B sites alternating with half- 
filled A sites is not possible. 

B. Reciprocal space charge ordering 

We now turn our attention to the charge populations 
of the one-electron (-) and (-1-) hybrid bands obtained 
in Eq. (O for J7 = 0. The upper (lower) tight-binding 
band is half-filled (filled) for any non-zero A in the t > 
caso^ while for t < this is only the case for: A > 
0.68|t|. Simple arguments might then suggest that if U 
is sufficiently large then the half-filled band may undergo 
a Mott insulator transition. However, for such a single 
band argument to be valid the half-filled + band must 
be sufficiently high in energy above the fiUed-band so 
that interband transitions induced by U can be safely 
neglected. 

In the basis of the non- interacting -|-,- band states the 
Hamiltonian is 



kQ.(T ^ k,k' ,qQi .Q9 ,a3 .Q4 

Q6) 



where the a's refer to the two values + and — and V is 
the Coulomb matrix describing the 16 different scattering 
processes between the bands: ec(k) of Eq. ([5]). The 
Coulomb matrix is 

y(kiai,k2a2,k3a3,k4a4) = (17) 
f//3(kiai)/3*(k2a2)/3(k3a3)/3*(k4a4), 



with 



/3(k±) = 



1 



a=Fk ^k± — ^k=F 



(18) 



where ^k± is given by (|4]) and ak, the normalization 
constant of Eq. ([3]). 

The occupation of the non-interacting bands 
is obtained through the expression: n± = 

Ek.(*o|cJ±,Ck±.|*o)- 

In order to investigate the populations of the non- 
interacting bands at large U wc plot Lanczos results 
for n_ - n+ for N,, = 18 in Fig. M For U = 2\t\, 
we find that n_ — = 1. However, 7i_ — < 1 
for large U. The interpretation of this result is compli- 
cated as n_ — ri+ conflates two effects: (i) charge transfer 
between the bands, effectively doping the 4— band with 



electrons from the lower band; and (ii) for U ^ the 
bands are no longer eigenstates, thus the physical inter- 
pretation of n_ — n+ is unclear for large U. In spite of 
this interpretative difficulties it is interesting to note that 
the behavior seen in Fig. [H] differs from a recent mean- 
field approact^ which includes local electron correlations 
only. However, at present, it is not possible to conclu- 
sively determine whether this is because non-local elec- 
tron correlations may play an important role or because 
of the strong interband scattering induced by the large 
U which will eventually destroy the reciprocal space de- 
scription. This question is particularly important given 
the proposed role of tiny hole densities in the -f-band^ in 
explaining the apparent discrepancy between the insu- 
lating behaviour suggested by resitivity^^i^ ARPES^i 
and optical conductivity^ and the observation of metal- 
lic quasiparticles via Shubnikov-de Haas experiments^ 
on Nao..5Co02. This question requires future investiga- 
tion. 
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FIG. 7: (Color online) Charge disproportionation, ns — wa, 
between inequivalent rows in the ionic Hubbard model [1] In 
the top panel tight-binding (dash-dotted lines) exact results 
for U = Q are compared with Lanczos diagonalization (open 
symbols) for the Ns — 18 tilted cluster of Fig. [3]showing good 
agreement. In the bottom panel the dependence of ns ~ ua 
with U is shown from Lanczos diagonalization for t < on 
the same cluster. 



C. Charge gap 

In order to understand the electronic properties of 
model (H)) we now discuss the the charge gap and its de- 
pendence with charge order driven by A and the Coulomb 
repulsion. The charge gap is defined by equation ([2]). 
When t is small, the lowest order correction to excita- 
tion energies come from the kinetic energy gain due to 
the propagation along the B(A)-chains of a hole (dou- 
blon) when extracting (adding) an electron to the zeroth 
order ground state configuration. Using degenerate per- 
turbation theory on Hamiltonian © the gap to 0(i^/A) 
is 



A, = A-2|<| 



8tVA - 2^^^ - 8tVA' 



(19) 



where SE^j/ is the energy change when adding a hole to 
a single half-filled A-chain modelled by the t — J model 
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FIG. 8: (Color online) Difference in filling of the hybrid ±- 
bands as a function of A/\t\ for t < and several values 
of (7/|t|. Results are from Lanczos calculations on 18-site 
clusters. 



with J given by Eq. ([7]). The second order hopping 



term t 



(2) 
f// 



2t^/A, contributes to the propagation of an 
antiholon in the A-chain. In the antiferromagnetic case, 
J > 0, the energy, &E\^^ in per A-chain site is given 
by the Bethe ansatz expression)^ 



-2\t\ 



sin(nyi7r) 



8J(ln2)n^ 



1 — sin(2n^7r) 
2nyi7r 



J In 2, (20) 



where riA = 1 — ^/N^ is the number of electrons in a 
single A-chain of iV^ sites when a single hole has been 
added to the otherwise half-filled chain. In the ferromag- 
netic case ( J < 0) 



El^' = -2\t\ - J/2. 



(21) 



Hence, the charge gap is found to be larger for t < 
than t > 0, as shown in Fig. [51 due to the geometrical 
frustration. In contrast, on the square lattice, Ac, does 
not depend on the sign of t. The dependence on the sign 
of t becomes even more apparent for A ~ \t\, where it is 
clear that for f < the gap is significantly larger than 
for i > 0. Indeed, it may be that for small A/|i| and 
t < the system is an insulator and for t > it is metal- 
lic. However, finite size effects prevent us from making 
a definitive statement about the existence of a metallic 
state for small A/|i|. We note that dependence on the 
sign of t is the opposite from what one would expect from 
weak-coupling arguments. For t > and U = A = and 
at 3/4 filling the Fermi surface has perfect nesting and 
there is a van Hove singularity in the density of states 
at the Fermi energy. Hence, weak-coupling arguments 
would suggest that for this sign of t the system would 
have a greater tendency to density wave instabilities and 
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FIG. 9: (Color online) Dependence of the charge gap on the 
on-site potential A. The charge gap for I] = 100|t| on a 
Na ~ 18 tilted cluster (top) and a A^s = 16 ladder-type cluster 
(bottom) (see Fig. |3]) is shown. Dashed lines denote the 
results of the strong coupling expansion (|19p for comparison. 



insulating states. Further, our mean-field slave boson 
calculations^ also predict that the insulating state is more 
stable for t > than i < 0, in contrast to the results re- 
ported in Fig. O 



D. Magnetic order 

Spin correlations in the model are analyzed through 
the static spin structure factor: 



^(q) = ]^ E exp^''-^'''"''^' < S^S] > (22) 



where = (n^i — njj^)/2 is the z-component of the spin 
at the lattice site Rj. 

The dependence on A of the spin structure factor, 
S'(q), is shown in Fig. [10] for two different wavevectors 
with U ~ 100|<|. For t < the results indicate a tran- 
sition to the magnetically ordered state with wavevector 
Qi = (tt, 7r/\/3) (which implies that the spins within the 




FIG. 10: (Color online) Development of magnetic order. The 
static spin structure factor S{q) at the wavevectors (a) Qi = 
(0, 7r/\/3), associated with C-type antiferromagnetisro^^ (see 
Fig. [T] left panel), and (b) Q2 ~ (7r,7r/\/3), associated with 
G-type antiferromagnetisro^ (see Fig. [1] right panel)), are 
shown for U — 100\t\ as a function of A, the difference be- 
tween the local energies on the A and B sublattices. Calcula- 
tions are performed on the Ns = 18 cluster. 



A-chains are antiferromagnetically ordered and antifer- 
romagnetically coupled with neighbouring A-chains) at 
about A = 2\t\. This is the spin pattern shown in Fig. [T] 
For t > there is a range of A for which the A spins are 
ferromagnetically coupled. The FM region obtained from 
Lanczos diagonalization is in agreement with the condi- 
tion J < extracted from a strong coupling expansion. 
This condition gives a FM region for 5i < A < ^/2U. 

The wavevector of the magnetic order and the associ- 
ated magnetic moment of the ground state of the model 
with t < and f/ » A > |t| (see Fig. [T]) are consis- 
tent with the antiferromagnetic ordering wavevector ob- 
served in Nao.sCoO ^^^i"^^ . However, the observed charge 
transfer between A and B chains in Nao.sCoO a^^'^^i'^^ is 
much smaller than the complete charge transfer sketched 
in Fig. [TJ Such a large magnetic moment in the pres- 
ence of a weak charge transfer between A and B sites is 
not expected from either classical or weak coupling argu- 
ments. Nevertheless, our previous exact diagonalization 
calculations for U > W (see Fig. 4 of Ref. 0) found 
a substantial magnetic moment and small charge trans- 
fer, even for small A, consistent with the experimental 
results. 
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FIG. 11: (Color online) Density of states of the charge transfer 
insulator. We take U = 100|t|, A = 10|t| and t < 0. The inset 
shows the low energy part where the charge gap agrees with 
the strong coupling expression Equation (|19|) . The excitation 
energies in the atomic limit {t = 0) are shown by the arrows 
below the abcissa. The chemical potential {lj = /i) is shown 
by the vertical dashed line. 



IV. DYNAMICAL PROPERTIES 

In this section we discuss dynamical properties of the 
model (HI). We use Lanczos diagonalization on Ng = 18 
clusters to compute the one-electron spectral density and 
the frequency dependent conductivity. 



A. One-electron spectral density 

The spectral density per spin is 



^ |(*„.(iV - l)|c,,|«'o)p5(a. + {Ern{N - 1) - Eo{N))) 

in 

+ |(vI/,„(7V + l)|c+ |vI/o)|2^(c^ - (£;,„(iV + 1) - E,,(N))), 

(23) 

where Em{N ± 1) is the spectra of excitations of the full 
quantum many-body problem with A'^ ± 1 electrons and 
\'^m{N zk 1)) its associated wavefunctions. Eq{N) is the 
ground state of the N electron system with wavefunction 
l*o>. 

In Fig. [TT]wc show the spectral density for U = 100|t| 
and A = 10|t|. As the system is well into the strong 
coupling regime, U » A >> \t\, we can understand 
the main excitations observed on the basis of the atomic 
limit (i = 0). 

For t = extracting an electron from the lattice can 
lead to three possible excitation energies: — 3A/2, A/2 
and —A/2 + U . While adding an electron we only have 
one excitation energy at A/2 -|- U . Therefore, four peaks 
at these excitation energies are expected in A{lj). A 




FIG. 12: (Color online) Evolution of the energy dependence 
of the density of states with decreasing charge transfer. We 
take U = 15|t| and vary A from 10|t| (top, charge transfer 
insulator) to A = \t\ (bottom, covalent insulator), all with 
t < 0. The chemical potential [uo = /i) is shown by the vertical 
dashed lines. 



gap of Ac = A is therefore obtained when t = {) typi- 
cal of a charge transfer insulator. When the hopping is 
turned on, the lowest order correction to excitation en- 
ergies comes from the propagation of a hole (doublon) 
along the B(A)-chains when extracting (adding) an elec- 
tron to the zeroth order ground state configuration. As 
the hole (electron) added can be on any site of the chains 
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the ground state is NA{NB)-io\d degenerate. This degen- 
eracy is Ufted at the first order in t, where the excitation 
energies for removing an electron are A/2— 2|t|(sin(fc) — 1) 
and A/2 — 2|t| cos(fc), from the A and B-sites respectively. 
Adding an electron to the A-chain leads to a doublon 
with excitation energy A/2 + U +2\t\{sm{k) — 1) Thus, a 
characteristic one-dimensional broadening of 4\t\ to the 
four peaks should be expected and the gap is reduced 
from the atomic limit result to Ac = A — 4|t|. The hy- 
bridization between the chains lowers the ground state 
energy for the N electron configuration, due to virtual 
excursions from a B-site to a nearest-neighbour A-site, 
by — 4|fp/A with no cost in Coulomb repulsion energy. 
Therefore, the lowest TV -I- 1 electron excitation energy 
is pushed upwards by -t-4|tp/A and the — 1 down- 
ward by — 4|ip/A. This leads to an increase of the gap: 
Ac = A - 4\t\ + 8f/A. The final charge gap, Ac, includ- 
ing the higher order corrections of Eq. coincides 
with the numerical calculation shown in the inset of Fig. 

m 

In Fig. [12 we show the evolution of A{uj) with A for 
U = 15\t\. The four peak structure discussed above for 
the CTI remains for this smaller value of U and A = 10|i| . 
As A decreases the peaks broaden due to hybridization 
between the A and B chains and shift in energy. For 
A = |t|, A{oj) contains a lower Hubbard band (LHB), 
an upper Hubbard band (UHB) and most of the spectral 
weight is around the chemical potential lu = ^. We find 
that the energy difference between the LHB and UHB is 
much larger than U , which we attribute to hybridization 
between the chains. Thus we identify this regime as a 
covalent insulator (CVjJ^ 



B. Frequency dependent conductivity 

The incoherent part of the optical conductivity is cal- 
culated through the current correlation function: 



0.1 



|(vI/,„(iV)|j.|vI/o)P 



Eri 



■d{co - (£;,„ - £^0)), 
(24) 

where, is the a;-component of the current operator, 
j = Z]<ij>,7,<T(R-i~R'i)cici+^o-, and R; is the position 

th 

of the i lattice site. All nearest- neighbor sites entering 
the sum in the current are denoted by 7. 

The evolution of the optical conductivity, cr(w), with 
A is shown in Fig. 1131 There are two main absorption 
bands in the optical spectra. One is fixed at large energies 
of about U and is associated with excitations between the 
Hubbard bands and a lower band which shifts with A. 
The lower absorption band is due to excitations associ- 
ated with transferring an electron from the — band to 
the -|- band. These produce a continuum of particle-hole 
excitations of width of order W ^ 9\t\. 

We have also calculated the Drude weight for a range 
of parameters in the model Hamiltonian. We do not show 
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FIG. 13: (Color online) Frequency dependence of the optical 
conductivity for several values of A/t and fixed U = 15\t\ 
with t <0. 



the results here because due to finite size effects the de- 
tailed interpretation is not clear. However, the trend is 
clear: as U/\t\ and A/\t\ increase, the Drude weight de- 
creases significantly. 



V. CONCLUSIONS 

We have considered the electronic properties of an ionic 
Hubbard model at 3/4-filling with stripes of alternating 
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on-site potential. This model has a rich phase diagram, 
in which various types of insulating and metallic states 
compete. A charge transfer-type insulator, a Mott in- 
sulator and a covalent insulator occur in different U-A-t 
parameter regimes. The geometrical frustration of the 
triangular lattice leads to different magnetic properties 
depending on the sign of t. For i < an antiferromag- 
netic interaction occurs whereas for t > a ferromagnetic 
coupling occurs in a broad range of parameters. 

At U >> A >> a charge transfer insulator of 
doubly occupied chains of i3-sites alternating with singly 
occupied A-chains occurs, i.e., ub — nA = 1- As A is 
decreased the system remains insulating although the 
charge disproportionation between sites is incomplete: 
Ub -riA < 1- 

The insulating state of Nao.5Co02 is characterized by 
a small charge modulation, a small charge gap, and 
strong Coulomb interaction. Electronic structure cal- 
culations suggest that Nao.5Co02 is in the parameter 
regime, U >> \t\ and A ~ \t\ which is different to the 
strong coupling parameter regime. This regime is diffi- 
cult to analyze within weak coupling perturbation the- 
ory and numerical approaches arc helpful. The model 
with no charge modulation, A = 0, reduces to a highly- 
doped Hubbard (or t — J) model on a triangular lat- 
tice which is believed to be metallic. However, under 
a weak external periodic potential, A ~ 0(|f|), our ex- 
act diagonalization analysis suggests a non-zero charge 
gap. Unlike the charge transfer insulators proposed pre- 
viously for Nao.5Co02 this insulator is characterized by 
small real space charge transfer: ub — ua « 1- An 
insulating state induced by strong hybridization of the 
non-interacting bands is realised and is reminiscent of 
the covalent insulator (CI) description of some transition 
metal oxides 4^ 

The behavior of the charge gap for U » \t\ has been 
further explored by an exact analysis of two and four- 
site clusters which indicates that the charge gap is al- 
ways enhanced with A even at small values. A different 
dependence on A of the ground state many-body ener- 
gies Eo{N), Eo{N ~ 1) and EniN + 1) is found which 
reflects the different nature of the bonds formed between 
inequivalent A and B sites. Two-electron bonds are 
well described as valence bonds which contain the effects 
of strong electronic correlations whereas three-electron 
bonds are accurately described by "molecular" orbitals 
which are uncorrclatcd. While the former type of bond 
depends weakly on A the latter docs not. 

Optical conductivity experiments suggest a gap of the 
order of 0.020 cV |t|/5), a sharp peak at about 0.026 
eV (~ |t|/4) which is at the lower edge of a continuum 
of excitations which reaches energies up to about 0.9 cV 
(~ 9|t|)^. This behavior is consistent with the low en- 
ergy adsorption band found in the calculated (t{uj) (see 
Fig. [T3)) which is located at about A ^ 0(|t|) with the 
continuum being the whole set of particle-hole excitations 
between the two hybridized bands which spread over the 
whole bandwidth ^ W. 



Recent experiments on Na2;Co02 for x — 2/3 the Na 
ions induce a charge ordering pattern with filled non- 
magnetic Co'^''" ions arranged in a triangular lattice and 
Qq3.44+ magnetic sites forming a kagome lattice structure 
with the transferred holes moving on it. These experi- 
ments are important as they relate the charge order with 
the different magnetic and electronic properties of the 
material. The present ionic Hubbard model modified to 
include such ordering patterns could be used to explore 
the unconventional metallic properties of Na2.Co02 at 
X = 2/3. 

We now briefly discuss the relationship between the 
results we obtained here and those we recently obtained 
for the same model with a slave boson mean-field theoryi^ 
Slave bosons give an insulator only for A > 8\t\ (A > 5|t| 
) for t < (i < 0) whereas exact diagonalisation suggests 
that the ground state is insulating even for A ^ 

An important open question that this study raises is, 
what is the ground state for small A/t? The tempera- 
ture dependence of the magnetic susceptibility of the t~J 
model on the triangular lattice has been calculated using 
exact diagonalisation on small clusters4^ The bottom left 
panel of Figure 6 in Ref. shows that for J = (i.e., 
L/ ^ oo in the Hubbard model) that at 3/4 filling that 
for all temperatures above about 0.4|i| that the suscep- 
tibility is the same as that for localised non-interacting 
spin- 1/2 particles. The susceptibility has a maximum at 
about 0.3\t\ and then decreases with decreasing temper- 
ature to a value about 2-3 times the value for U ~ 0. 
These results raise the question as to the nature of the 
ground state and the tendency of the electrons to be- 
come localised and the spins to antifcrromagnctically or- 
der, even in the absence of an exchange interaction, due 
to kinetic antiferromagnetism4^ 



APPENDIX A: HEISENBERG EXCHANGE 
COUPLINGS 

In this section we discuss the various contributions to 
the nearest neighbour exchange coupling, J, between the 
A-sites. Taking \'^o > as the ground state configuration 
in the strong coupling limit U » A » t, there is no 
correction to the lowest order in the kinetic energy. To 
0{t^) we have the usual superexchange antiferromagnetic 
contribution: J = AP/U. To 0(<'^) "ring" exchange pro- 
cesses around the 3-site plaquette of the type shown in 
Fig. [T?] remove the spin degeneracy and were already 
discussed by Penc and collaborators [i^ in a Hubbard 
model on a zig-zag ladder. The energy of the singlet state 
in Fig. [14] is shifted by 4t'^/A^ while the triplet state by 
—At^/A'^. These shifts are opposite to the two site case. 
As there are two possible ways of going around the trian- 
gle in Fig |14l and there are two neighbouring B-sites (one 
below the two A-sites as shown in Fig. [14] and another 
above) the final contribution to the effective J at 0{t^) 
\s J ~ Et — Es ~ — 8i"^/A^ enhancing the ferromagnetic 
tendencies as compared to the ladder caso^ by a factor 
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FIG. 14: (Color online) Three site "ring" exchange processes 
contributing to the exchange interaction, J, between neigh- 
bour spins in an A chain of the t — J — J± model (compare 
Equation Q). 




FIG. 15: (Color online) Four site "ring" exchange processes 
of 0((t/A)'* contributing to the exchange interaction between 
two neighbouring sites in the A-chains in model ([6]). 



all possible exchange processes, a total contibution to J 
valid to ©((t/A)*) is 

4t2 8t^ 16t^ 48*4 IGf^ 



APPENDIX B: GROUND STATE 
WAVEFUNCTIONS FOR THE FOUR-SITE 
CLUSTER AND VALENCE BOND THEORY 

In this appendix, we discuss exact ground state wave- 
functions on the Ng = 4: cluster of Fig. [4] with t' = t 
and = 5, 6 and 7 electrons. Valence bond (VB) states, 
which are neutral configurations formed by two neigh- 
bour electrons in a singlet are found to describe the 
exact wavefunction accurately for U >> \t\ and t < for 
N = 6. This is because ionic configurations have a neg- 
ligible weight in the full wavefunction at large-C/ values. 



1. Ground state for N = 6 electrons 



of two. In contrast "ring" exchange processes around a The exact ground state wavefunction, |^'o(6) >, is well 
4-site plaquette {0{t'^)) of the type shown in Fig. [15] described by the resonance between different possible VB 
lead to an AF contribution to J = 40t"'/A'^. Including states between electrons inside the cluster: 



|«'o(6) >«a 








Ti 

I n 

T 



(.Bl) 



Here, and below, VB states are colored red. Note that the horizontal sites are A sites while the vertical sites are B 
sites as sketched in Fig. [H For f7 = oo and A = 0, the energy of the RVB state is E^^^{6) = -3.3723|t| and the 
weights of the wavefunction are a = b = 0.4544, and c = 0.3831. This RVB wavefunction gives an accurate description 
of the exact ground state energy which for U = 100|t| is: Eq{6) — — 3.453|<| and a wavefunction described by (|Bip 
with a ~ 0.459, b = 0.445, c ~ 0.382, plus small ionic terms. Note that, on the four site cluster, the A and B sites are 
not equivalent even for A due to geometry of the cluster (cf. Fig. 0]) 

The RVB wavefunction also accurately describes the ground state energy in the limit: U >> A >> \t\. The energy 
of the RVB state for A = 10|t| is E^^^ = -10.433|i|, and the wavefunction has weights a = 0.977, b = 0.021, and 
c = 0.1058. This is in good agreement with the exact ground state, which has Eo{6) = — 10.453|t| and a = 0.976, 
b = 0.024 and c — 0.108 (plus small ionic terms). Thus the VB formed between the two A sites dominates the 
wavefunction. 



2. Ground state for N — 7 electrons 



The 7-electron system contains one hole which can hop around the cluster so the Coulomb interaction has no effect. 
The ground state is simply a linear combination of the states with one hole in the cluster: 
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(B2) 



Note that the e and d coefficients are different even for A = due to the geometry of the cluster, cf. Fig. SI 
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All many-body configurations contain VB singlets ex- 
cept for the last term (proportional to h). The electrons 
in the A and B sites prefer to align with antiferromag- 
netic order in contrast to predictions of RVB. This is due 
to the different number of spin up. 2, and down elec- 
trons, 3, for TV = 5. When a spin up is located at an 
A site, the remaining two spin down electrons at the B 
sites gain energy by aligining AF with the A electron. 
This is more favorable energetically than having a VB 
betwen two electrons at A and B sites and a down spin 
at the left B-site. For U = 100|i| and A = 0, / = 0.423, 
g = 0.326 and h = 0.462 with the ground state energy, 
Eq(5) = 97.590|t|. For A = 10|t| the wavcfunction is ap- 
proximately given by / = 0.691, g = 0.076 and h — 0.107 
with ground state energy, Eo{5) = 94.615|i|. 



.(B3) 
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3. Ground state for N = 5 electrons 

The 5-electron ground state wave function is well described by the resonance between the following states: 
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